Detecting stable massive neutral particles through particle lensing 
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Stable massive neutral particles emitted by astrophysical sources undergo deflection under the 



I. INTRODUCTION 



?—( , gravitational potential of our own galaxy. The deflection angle depends on the particle velocity and 

^2 ' therefore non-relativistic particles will be deflected more than relativistic ones. If these particles can 

be detected through neutrino telescopes, cosmic ray detectors or directional dark matter detectors, 
their arrival directions would appear aligned on the sky along the source-lens direction. On top 
Vi ■ of this deflection, the arrival direction of non-relativistic particles is displaced with respect to the 

CN ' relativistic counterpart also due to the relative motion of the source with respect to the observer; 

this induces an alignment of detections along the sky projection of the source trajectory. The final 
alignment will be given by a combination of the directions induced by lensing and source proper 
motion. We derive the deflection- velocity relation for the Milky Way halo and suggest that searching 
V^ ' for alignments on detection maps of particle telescopes could be a way to find new particles or new 

—7 , astrophysical phenomena. 

O 

C^ ' As it is well-known (see e.g. [16]), a light ray trajectory is deviated by a compact mass M by an 

amount 

>: a = '-f (1) 

^: 

T- 1- ' where Rs — 2GM/c^ is the Schwarzschild radius of the deflector (or "lens") and £, is the impact 

l/~j , parameter, i.e. the shortest distance of the particle to the deflector along the imperturbed path. 

The angle a, denoted lensing angle, defined as in Fig. ((T]), lies in the plane observer-lens-source. 

}^ [ This equation holds true for small lensing angles and for a deflector size much smaller than the 

^SJ ' source distance. 

The analogous deflection for a non-relativistic (NR) particle is easy to derive in the Newtonian 
limit. If a particle is emitted by an astrophysical source with a non-relativistic velocity v (in units 
of c), its trajectory will be deflected by any mass close to its trajectory by an angular amount equal 

X : to (e.g. M P- 2) 

* = 4 (2) 

again in the limit of small deflection angles and deflector size. Crucially, this deflection is amplified 
with respect to the relativistic one by the v~'^ factor. Because of the strict analogy, we refer to this 
deflection as particle lensing, although in neither case there is any focusing of images. 

If the particle is neutral, its trajectory will not be further deflected by the galactic magnetic 
field and by estimating a, ^ and Rg we can infer the particle velocity and, given the kinetic energy, 
its mass. If particles of the same kind are emitted by an astrophysical source with a range of 
velocities, they will all travel on the same plane observer- lens-source and reach the observer with 
different angles a{v). The observer will therefore detect particles with different energies aligned 
on the sky in the direction source- lens (see Figs. I2l3p . In this paper we suggest that this peculiar 



alignment of detections in particle telescopes, as e.g. neutrino telescopes (e.g. jlSl) cosmic ray 
detectors (e.g. (l^l) or future directional dark matter detectors (e.g. [2,113)) could be used as a 
test to reveal the existence of stable massive neutral particles (SMaNPs) emitted by astrophysical 
sources inside or outside our Galaxy. Let us remark from the outset that although the method we 
propose is extremely simple, it is not easy to identify possible targets for its application within the 
current frontiers of particle physics and astrophysics and it is therefore to be meant as a search for 
unexpected phenomena. For this reason, we will not investigate in any detail how and where such 
neutral particles could be produced and accelerated. 

Since dark matter particles are also expected to be SMaNPs, it is clear that a novel way of 
detecting them would be potentially very interesting. However, the particle lensing method ap- 
plies only to particles emitted by astrophysical bodies and clearly cannot be employed to detect 
halo dark matter particles. Possible candidates are massive sterile neutrinos (e.g. [8]) or exotic 
components of cosmic rays (e.g. heavy axion-like-particles [9|, Q-balls [lO|0i neutral strangelets 
|12| etc.) emitted by high-energy sources like supernovae remnants [a| or microquasars [20|. In 
principle, also transient phenomena like supernovae or binary star collapse could be interesting 
sources; however non-relativistic particles will arrive scattered for thousands or million years after 
their optical counterpart. Only if the flux is sufficiently large, could one still detect alignments of 
slow neutral particles with just slightly different velocities with arrivals spaced by days or months. 
In this paper we do not consider transient phenomena and time delays and assume that a sufficient 
and continuous flux of particles reaches the Earth. 

Two bodies could represent useful lenses: the Sun and the galactic halo. As we will show below, 
the Sun in quadrature with respect to the source will deflect SMaNPs above 0.1 degrees if they 
are slower than 1500 km/sec (in the observer rest frame). Our galactic halo can deflect particles of 
the same velocity or smaller up to several degrees and its therefore the primary choice for the lens. 
Here "deflection" means in fact "deflection with respect to the optical or relativistic counterpart" 
of the source, or "differential deflection" for short. The angle of deflection a itself is in fact not 
observable since in general we do not know the unlensed source position. If the galactic halo acts 
as a lens, and assuming radial symmetry of the halo potential, the detection pattern would show 
up as straight alignments on the sky pointing toward the galactic center, as illustrated qualitatively 
in Fig. ([3]). Extragalactic sources and lenses are also possible but we will not consider them in this 
paper. 

There are three obvious problems with this technique. The angular resolution for non-relativistic 
particles is currently very poor: directional dark matter detectors exploit the fact that the deviation 
between the recoil direction of the detector's nuclei and the arrival direction is related to the recoil 
energy but for non-relativistic par ticles the resolution of current and planned detectors is very low, 
no better than 15 degrees [2, llTj . Cosmic ray experiments and neutrino detectors reach much 
higher resolutions, down to 0.1 degrees, but only for ultra-relativistic particles. In this paper 
we assume optimistically that the resolution of 1 degree can be reached in the near future also 
for non-relativistic scattering. The lower the resolution, the more difficult is to distinguish the 
lensing alignment from chance alignment of background sources, halo dark matter or other unlensed 
particles. 

The second problem is that the lensing technique is sensitive to an extremely narrow energy band. 
Typically, if the energy is larger than the particle mass by more than 10^^ or so, the deflection by 
the galactic halo reduces to the sub-degree level and cannot be resolved by current or foreseeable 
particle telescopes. In fact, since they must be non-relativistic, all detectable lensed particles of the 
same species should have almost exactly the same energy, practically coinciding with their mass. 
This reduces drastically the flux of non-relativistically lensed particles but on the other hand could 
simplify their search, since only equal-energy detections could have been aligned due to particle 
lensing (this of course applies only if the particles transfer all, or a fixed fraction, of their energy to 
the detector). 

The third problem is that, if the source is moving with respect to the observer, the non-relativistic 
particles that reach today the Earth were emitted when the source was in another position. This 
induces a spread of the arrival directions with respect to the relativistic counterpart. The particles 
will then appear aligned along the sky projection of the source trajectory ("proper motion vector"). 
This "proper motion alignment" (PMA) adds to the lensing alignment and in general produces a 




Figure 1: The geometric configuration of tlie lialo lens system. Tlie observer is in O, tlie source in S, its 
image in S', and the lens is the filled circle. Point C is added to ease the discussion. In the following we 
always assume a, a to be small. The angles 6, 6 are also small in Sect. IIIII while are left free in Sect. IIVI 



new alignment intermediate between the source-lens one and the source proper motion vector. If 
the source proper motion is known, the PMA can be disentangled from the lensing effect. This 
paper is devoted to the calculation of the lensing alignment; we confine in an Appendix a brief but 
sufficient treatment of the PMA. In general, the proper motion does not spoil the lensing signal 
and in principle can be accounted for. 

In this paper we derive the geodesic equation in a general-relativistic setting, bridging the gap 
between Eq. ([ij and © and we derive the angle of deflection as a function of velocity or of energy 
per mass. Although as we already emphasized we cannot at the moment point to any specific 
realistic target, the method we propose is so straightforward to carry out on existing and future 
detection maps that we feel it is an interesting addition to the astroparticle physics toolbox. 



II. GEODESIC EQUATION 

We begin by writing down the equations of motion in the weak and slow- varying field limit but 
without limitations on velocities. We assume the usual longitudinal metric for scalar perturbations 
in a Friedmann-Robertson- Walker 



ds^ = -(1 + 2*)rft2 + a^{l + 2<^)dx'dxi 



(3) 



The passage to a Minkowsky metric, appropriate for Galaxy lensing, will be carried out at the end. 
Defining the affine parameter r such that dr^ = —ds^ we can write 

dr ^ dt{l + *)r (4) 

where for small $, ^, 

r= [l-i;2(l + 2($-*))]i/2 

and where we define the peculiar velocity 

"" " "^ (^) 

We need to derive the geodesic equations for a particle 

^ + r>".^ = o (6) 



where the four velocity is 

Here we assume that the gravitational potentials are small, but the velocity can be relativistic. We 
assume from now on that the time derivatives of ^, $ are negligible, i.e. that the lens is a static 
object. Notice that d/dt is a total derivative and therefore 

divUl - *)/r) 1 , , , vv v'vi'i' , - v'^v^v^^ ,- 

dt r^ r2 ar2 > ^ > 

where w = |v| . 

From the geodesic equation with fi = i we obtain to first order in <&, $: 

(9) 

with 

e„ = 2*(1 - 27172) (10) 

and where we used the usual special-relativistic factor 

^' = {l-vT' (11) 

We put now $ = — ^, valid for ordinary matter. Then we obtain 

v^{l + ei) + J^v^vv{l + £2) = -Hv,{l + £i) + (272 + l)w,(i;, Vj*) - (1 + w^)^^^ (12) 

and from now on all spatial derivatives are with respect to physical distances r^ = axi. The fi = 
equation is 

^'^vv{l + £2) = -v'^H{l + £1 - 4*) + (27^ - 3)(wj Vj*) (13) 

(this can also be obtained up to terms higher order in ^ by multiplying Eq. (fT^ by w') and we can 
use it to simplify Eq. (fT2|) obtaining 

72v(l + £1) = -iJv(l - 2*) + 472v(v ■ V*) - (1 + 2w2^2)y^ ^^^^ 



However all the equations here assume that 7^^ -C 1, otherwise we would need to include terms 
like 7^\1/^ etc.. In this limit we can approximate ei « so 

72V = -Hw + 72[2vv • V(* - $) - w2V(* - $)] - V'J' (15) 

In Appendix A we generalize this equation to a scalar-tensor gravity. Although the limit v ^ 1 
cannot be naively taken from this equation, it turns out that it approximates the correct result. In 
fact, if a light ray passes near a spherical potential $, in the vicinity of the potential (where the 
effect is largest), the direction is orthogonal to the gradient and therefore v • V(^ — $)—)• and for 

7 ^^ 00 

V = -t;2V(* - $) = -V(* - $) (16) 

which is indeed the lensing equation (see e.g. (7|). 



III. PARTICLE LENSING 

We evaluate now the lensing of relativistic and non-relativistic particles. From here on, we restrict 
to a Minkowky metric, i.e. we put H = and as before we put $ = — ^. Eq. (jTS]) becomes 

72v = -272[2w V$-i;V$] + V$ (17) 

We integrate this equation along a slightly distorted path. We can assume that the real trajectory 
is only weakly distorted with respect to the imperturbed one, which is given by 

vdt — vadrj — adr (18) 

with V = const. Then on the rhs of Eq. (IT71) we can assume that v follows the imperturbed path. 
We consider a point-like mass (lens) located at a distance r^ from the observer (see Fig. [IJ. We 
have then the standard Newtonian potential 

We use a frame such that the particle is described by the coordinates {rOi , r02 , t) with small 
deviation angles Oi^2- On this path, the term v • V$ changes sign when the particle overcomes the 
lens and therefore its total effect for a spherically symmetric lens is zero. We are then left with 

(z = l,2) 

._^2d(^^2z;2v.$ + 7-2v,ci> (20) 

or 

^ = 2V.$ + 4^V,$ (21) 

dv^ u^7^ 

where now spatial derivatives are with respect to comoving coordinates. The total lensing can be 
obtained by summing the two terms on the rhs. Let us analyse them separately. 
The first gives the standard weak-lensing equation. Its solution is (see eg Q) 

a^^0^-0l^2f^^d/-^'^^^V,.^, (22) 

where Vj^ is the gradient along the angular directions 61^2 , ^s is the source distance and 01 is 
defined as 

- -^'^'-;^l ,^,JZ;Z,„.A n'y> ,24) 



where R' — {r'9[, r'9'2) and R = {tlSi, ^£^2) and where if the lens potential is sharply peaked at 
the lens position we can assume r = r^ in the slowly varying factors inside the integral (in other 
words we assume here that the lens is a point mass separated by the source by a small angle, i.e. 
the typical strong-lensing astronomical configuration; later on we generalize to a distributed mass). 
Finally we have (up to a constant independent of and therefore irrelevant when the gradient V_l 
is applied) 

(/.i(^)=4G^^^^^ /"d2^'I](R')ln|R'-R| (25) 

rsrL J 

where the surface mass density is 

S = /" drp{R, r) (26) 

The second term in (l?T|) can be treated similarly and gives 

M&) = 2G '/ 7 ""^ / d^R'j:{R') In |R' - R| (27) 

v^^-^rgrL J 

The second term is therefore (2w^7^)~^ times the first one, and dominates for w <C 1. The final 
result is 



a = V_L0 (28) 



where 



So finally we have 



2^2^2 TsTl 



4G(1 + TJITt)^:^ / '^'fi'S(R') In |R' - R| (29) 



1 ^ rs-rL f ,2^, S(^0 
2«V^ r, y'^ll-^' 



"(0 - 4G(1 + :^)^V^ / d'^ji^i^ C) (30) 



where ^ is a vector on the lens plane (impact vector). For a point lens of mass M, this gives 



2^2^2 ^^^ 



a = 4GAf(f + ;^-^^)^-^ (31) 



where ^ = |^| = r^ sin 5 is the impact radius. The angle a is related to the angle a (see Fig. [T|) 
between initial and final propagation (more often employed by astronomers) by the relation 

a = a^^^ = 4GA/(1 + -1^)1 = ^ (32) 

where Rs = 2GM/c? is the Schwarzschild radius and 

Eq. (I21D generalizes Eqs. (fTT^ . 

Using these formulae, we can employ the standard lensing equations, simply replacing T\Rs for 
Rs- Making reference to Fig. ([T]), for small (5, the imperturbed angle between source and lens, the 
ang ular separation between the direction of the lens and the direction of particle arrival is (see eg. 

m 

^i.2-^(^±v/4t2+,52) (34) 



where 

r = ni/2 S/l^l£ = Ti'/^ao (35) 

V rLTs 

and where ao is the standard characteristic angle. Finahy, the differential deflection, i.e. the 
angle between the different arrival directions of particles of velocity v with respect to the opti- 
cal/relativistic counterpart will be 

Ae = ^{yjmal + 52 - ^4al + S^) (36) 

where we only considered the solution with > 5 . From this equation we can identify three regimes 
of velocities. For very small w, the dominating term is 4naQ and we expect a 1/v behavior: 



M « nV2«„ .. ^-Jrs'-^^^ -. ^-J^ (37) 

V V rLTs vyrL 

(where in the last expression, here and in the next equation, we assume Vs ^ r^). For larger v, but 
still w ^ 1, we expect a l/w^ trend since 

Finally, for w — !> 1 , A0 quickly vanishes. 



IV. LENSING IN A HALO 

In the following we derive the lensing in our own galaxy halo, where the approximation of point 
mass and of small 8 is no longer acceptable. We assume a general spherical mass distribution M{r) 
centered on r^ = (r^ cos 6*, r^ sin 6*) on the propagation plane. In Fig. (HJ we display the geometric 
configuration of the problem. In this case the lensing equation amounts to the statement that 
S"S' + SC = S'C , from which, assuming that OA«rL cosi5 and OS'«OS, one finds to first order in 
a 

T — ' T T COS 

iaiiO ^ ia.n5 + a— (39) 

where if $ is the gravitational potential of the halo (see e.g. Eq. (US])) 

a{5,v) = 2n ^^ f' d/-^^^^^{\v~VL\)U=o (40) 

fs - tlcosO Jo Tgr dr] 



2n-^ 



tl sm f 



If" rs -r TLSmS d 

2 70 r^-r-LCOS^ rAI drj 



(41) 



where Rg — 2GM and M — A/ (r^ sin (5) and the distance between the lens center and a generic 
point r = (r cos 77, r sin ry) is 



■tl 



[r^ +rl- 2rrL cos{f] - S)] ^/^ (42) 



Notice that Rs is the Schwarzschild radius of the mass contained inside the impact radius. Notice 
also that here we put the coordinates of the center as r^ = (r^ sin 6, r^ cos S) , since at this order 
can be replaced by S. To the same order, we can approximate the particle trajectory with the 
horizontal straight line ry = (i.e. y = ). The dimensionlcss term inside the square brackets (let 
us call it the halo integral) 



H=- dr i ^ . — $ r - ri ^=0 

^ Jo i^s - TL COS 5 rM dr] 



(43) 



reduces to unity for A/ = const and for small 6. We have intact in this case $(x) = M /x and 

1 f^" . Tfi — r r^sinJ d f 1 



H 



2.0 
1 



dr- 



dr 







^ -r 
tl cos (5 



dr] V |r - ri, 



\n=o 



Ts ~ tl cos (5/ (r^ + r\~ 2rrL cos 5)^/^ 



(44) 
(45) 



The integrand is sharply peaked at r w r^ cos 5 and therefore we can replace r with r^ cos S in the 
first parentheses of Eq. (|45p and find, in the limit of r^ ^ r^, 



dr 



(rL sin (5)^ 



1 + cos S 



{r^ + rj^ — "irrL cos 6)-^^^ 



1 



(46) 



(the condition Ts 3> tl can in fact be relaxed; for instance, H approximates 1 even in the isosceles 
configuration r^ ~ 2rL cos S) . 

Since we have now a = 2IiHRs/{rL s\n9) we can solve the lens equation 



tan 6 = tan S 



sin 9 



where 



2Rs{rs -rLCOsS) 



TsTL 



(47) 



(48) 



is a minor generalization of the standard characteristic angle. In Eq. (|47p the factor 11 takes into 
account the velocity of the particles, the factor H takes into account the deviation from the point 
mass configuration. Both reduce to unity for the standard point-mass relativistic lensing. Although 
Eq. (|47p could be easily solved numerically, here we assume simply that cos 9 k, \ even when x is 
not very small, and obtain the solutions 



ii^2 = arcsm 



-(tan(5 ± V4t2 + tan^ 5) 



where now 



nuai 



Finally, the halo differential deflection is 

-(tan,5+y^ 



IS.9 = arcsin 



UHHal 



tan^ 6) 



arcsin 



- (tan 6 



AHal 



iai? 5) 



(49) 



(50) 



(51) 



which generalizes Eq. ((36)) . Notice that for tan (5 > 1 the argument of the arcsin exceeds unity; in 
this case one should solve the full Eq. (|T7)) . li H = 1 and moreover 5 <^ \ (along with ciq ^ 1, 
an approximation needed in both cases) , we recover Eq. p6p . The same three regimes of velocities 
discussed before can be found here, although with more complicated behaviors. 



V. SUN AND STAR LENSING 



For a star of 1 solar mass with r^ ^ r^ one has 



TL \ fL J 



10-5 f 1!11^ ] arcsec^ 



For the Sun at r^ — lAU we have then, using Eq. 



10^ 1 
At^s„„ = jarcsec 

^sun ^ 



(52) 



(53) 




Figure 2: Schematic representation of the deflection due to the galactic halo. 



so that for a source at S — 1 rad the deflection is 



A0,, 



10- 



4y2 



-arcsec 



(54) 



We will see in the next section that the solar deflection, although non negligible, is quite smaller 
than the halo deflection. Moreover, it can be evaluated and subtracted from the total deflection 
since we know position and mass of the Sun with great accuracy. 

Assuming a density of 10~^ stars per parsec'^, in a cone of 1 arcsec aperture and 10 kpc depth, 
one expect roughly one star. The deflection for particles passing 1 arcsec away from a 1 solar mass 
star lens at 1 kpc is 



^f^star — 



IQ- 



1 



10- 



2^272 



r arcsec 



2w2 



-arcsec 



(55) 



again using using Eq. ([55)) . The deflection from single stars is therefore negligible with respect to 
the solar one. To estimate the resulting effect of many stars we can integrate the random deflection 
variance (A^^tar)^ over a cone of aperture Omax (expressed in arcseconds) along whose axis the 
particle travels from the source to us. Then we have 



{^Ototf - 



27r 

C2 



dr 



de 



RlC' 



UsT 



(56) 



where C 
have 



2 • 10^ converts from radians to arcsec and Us — O.lpc ^ is the star density. Then we 



(A^t, 



27r rg 
C2"^3 






log- 



3-10 



-6 



arcsec 



(57) 



where in the last step we assume, for instance, dmax = 10 deg and 6„iin — 1 arcsec. The total 
random deflection induced by the stars along the line-of-sight is therefore negligible as well. 
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Figure 3: A depiction of the detection map witli non-relativistic lensing for a source witli negligible proper 
motion, purely for illustrative purposes. The little red dots represent the arrival directions of massive 
particles emitted by source 1 and 2. 



VI. GALAXY HALO DEFLECTION 



Let us consider now the amount of deflection induced by the Milky Way halo potential. We 
assume that a source lies at distance r^ from us in a direction that makes an angle S with respect to 
the Galactic center (GC). The GC is then at coordinates {Re cos 5, Re sin S) on the plane observer- 
GC-source if Re is the GC distance from us. 

As a first approximation we can simply take Eq. ([55]) and put M ~ 10^'^ Mq. In this case the 
characteristic angle is of the order of 



10^ 



Ikpc 



M 



(58) 



For w <IC 1 we have from Eq. (|38 



A0„ 



an 



1 



5 2z;2^2 



10^ 



, \kpc. 
■( jarcsec 



TL 



(59) 



Assuming 5 = 0.25 rad and r^ = 10 kpc this is w O.l/u^ arcsec. As anticipated, this is one order 
of magnitude larger than Sun's deflection in Eq. (j54p . This is larger than 1 deg for v < 0.005. For 
V ^ 0.001 the deflection is of the order of tens of degrees and the small-deflection approximation 
adopted throughout would fail. 
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Parameter 


Halo 


Sun 


Star 


M 


Mnfw 


Mq 


Me 


rL 


Re 


1 A.U. 


1 kpc 


Ta 


2rLCOs S 


10 kpc 


>rL 



Table I: Parameters used for the Milky Way NFW halo, sun and a typical Milky Way star. The distance 
of the Galactic Center from Earth is Re — 8.5 kpc. Mmfw{x) is calculated from the integral H43p as a 
function of the distance r from the center of the lens. Alternatively, for the Galaxy Halo, we also compare 
with a constant mass lens of mass M — MifLsinS). 



For a more precise calculation, we need to integrate Eq. (PT|) with r^ — {Re cos S, Re sin 5). The 
mass M(x) inside a distance x from the center can be taken by integratin g a model for the Galaxy. 
Here, just to produce a rough estimate, we employ a Navarro-Frenk- White [1J| profile for the galatic 
halo 



PNFw{r) 



xIpo 



r{Xs + ry 



(60) 



with parameters Xs = Vy/c where ry is the virial radius and c is the concentration. For our source 
we choose c = 12.5, r„ — 200 kpc corresponding to Xg — 16 kpc. The density parameter is 
po = 0.014Mq/pc^ is chosen such that the mass within the virial radius is AI — lO^^M© |3l. 1 18l| . 
Then we obtain 



MNFwix) =47r / y PNFw{y)dy 



and the potential 



$ 



NFW 



(x) 



inpoxllogil + x/xs) 



(61) 



(62) 



Then we can evaluate the halo integral H and use Eq. ((5T|) to calculate the differential deflection 
A6. The results are plotted in Figs. (|H) and ([S]). As can be seen, the halo lensing can be larger 
than 1 degree for v < 0.004. The dependence on 6 in the range 0.1-0.5 rad is quite modest, since 
for small d both the impact radius and the halo mass within the particle trajectory decrease. In 
Fig. ([6]) we plot the deflection versus E/m « 1 + w^/2. This shows that as soon as the particle 
energy is 10~^ larger than the particle mass, the deflection decreases to less than 1 deg, making it 
impossible to distinguish from the relativistic counterpart. The particle lensing method is therefore 
sensitive only to a very narrow energy band, AE/E w 10^^. 

The assumption of a smooth NFW profile is of course very rough. A realistic galaxy model should 
include at a minimum the disk and the bulge. As long as they are both radially symmetric, the 
particles will still be aligned along the GC direction (for the disk cilindrical symmetry is sufficient, 
provided also the lens belongs to the disk). Deviation from symmetry will induce deviations from 
the straight alignment. 



VII. CONCLUSIONS 



In the last few years we have seen the beginning of particle astronomy, i.e. the detection of 
particles of astrophysical origin with detectors that are able to reconstruct the arrival direction. 
Present neutrino and cosmic ray detectors already reach an angular resolution for high-energy 
events of 1 degree and better [1^ [l^; although the directional dark matter detectors that are 
currently investigated only reach down to 15 degrees or so [3, one can reasonably expect that 
future technology will sensibly improve this, especially if additional science cases can be found in 
support for it. 
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Figure 4: Differential deflection A9 in units of arcsec as a function of velocity, in units of the speed of light, 
for Milky Way halo, Sun and star. We show the A9 corresponding to 5 = 0.5 rad for a NFW potential 
4i{x) (solid red) and for S = 0.1 (solid blue). For S — 0.5 we plot also the curve assuming a constant mass 
M* — M{rLsinS) (dashed red). The curve for S = 0.1 shows clearly the three regimes v~^ ,v~^ and 7~^, 
from small to large velocities. The Sun (solid orange) is shown for S = 0.5 rad. A typical star in the Milky 
Way of mass M = Mq, located at vl = Ikpc and r^ ^ r^ with 5—1 arcsec is also shown (dashed orange). 
Resolution thresholds of 1 deg and 0.1 deg are overplotted for reference (solid black). 
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Figure 5: Zoom of Fig. (|4]). 
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Figure 6: Differential deflection versus energy/mass for the NFW halo with 5 = 0.1 and 0.5, same parameters 
as in Fig. ^. 

In this paper we explore the possibility of detecting stable, massive, neutral particles (SMaNPs) 
by using particle telescopes. We point out that if these particles are non-relativistic they will appear 
to an observer aligned in the sky due to the combined effect of non-relativistic lensing and of the 
source proper motion. This alignment is a very distinctive signature that seems difficult to be 
mimicked by background noise. The method can be applied to particles of any mass, provided of 
course their kinetic energy is above the detector's threshold. Searching for alignments in detection 
maps could then be a simple and effective way to search for new physics and new astrophysics. 

We derived the deflection-velocity relation due to halo lensing and to proper motion and found 
that particles with velocity around 1000 km/sec are the best target for detection. We also pointed 
out however that the energy bandwidth is extremely small, roughly AE/ E « 10~^ and therefore 
the useful flux from any given source is very much reduced. There are several extensions of the 
standard model that predict the existence of SMaNPs [^, [sl-llOl. [l3 | ; how and where these particles 
could be emitted by astrophysical bodies is of course a difficult problem that we leave to future 
work. 
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Appendix A. Geodesies in a scalar-tensor gravity 

For generality, we can include the effect of a fifth-force mediated by a scalar field, as in a modified 
gravity scenario (see e.g. |^]), by confornially transforming the geodesic equation. If we put 

9^.u = e^f^^g^, (63) 

we obtain 

Kx = r^A - m^Sii + ^^uS'i - ^'^g.x) (64) 

and 

uf" = e'^'^u'^ (65) 

Then, by assuming also that 93 = (50 is taken at first order and ip is negligible we have 

w{l+ei)+-i'^wvv{l+e2) = -i/"v(l+ei)-V(*-/3(/j)-u2V(*+^.^)+v(v-V(*+/3(p)) + 27^v(v-V*) 

(66) 

where H = H{1 — j3fj) . This reduces to the standard case 

V = -Hw - V(* - I3lp) (67) 
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for small velocities. The /i = equation is 

72^(1 + £2) = -Hv^{l + ei- 4*) + (v • V((272 - 3)* + ^ip)) (68) 

which again can be used to simplify Eq. ([55)) : 

-/H{1 + £1) = -Hv(l - 2*) + 4-/^-v{v ■ V*) - 2j'^v'^V^ - V(* - (3ip) (69) 

If we keep the two potentials distinct, we obtain 

j^<r{l + £1) = -iJv(l - 2*) + 72[2vv • V(* - $) - v^V(* - $)] - V(* - /3ip) (70) 

where 

£1 = 2^(7^ - 1) - 2*7^ (71) 

It is interesting to note that the motion depends on ^ — j3(p alone when slow, but on ^ = $ — vjf 
(and not on (p, due to the conformal invariance of electromagnetism) when ultrarelativistic. 

Appendix B. Flux from supernovae remnants 

We can estimate the possible order-of-magnitude flux of SMaNPs from a typical supernovae 
remnant (SNR) in the following way. The 7-ray energy flux from typical SNRs, e.g. the sources 
RX J0852.0-4622 and RX J1713. 7-3946, as observed by the H.E.S.S. Cerenkov detector P, [IH in 
an energy band AE is 

Nae = ^AS (72) 

ah, 

where 

g.lO-"™-,-.re,.-(^)" (73) 

is the differential flux in the range lOOGeV — lOTeV and F « 2. Both sources are estimated to lie 
at a distance 1 kpc, with a large uncertainty. The flux from similar SNRs at distance D could be 
taken therefore to be 

Let us now suppose that the flux of SMaNPs is a fraction e of the flux in 7-rays. In an energy band 
AE/E = 10^"* (which as we have seen is the typical bandwidth for particle lensing) we have a flux 

M^,-.W-'.,n-'yr-'{^y'*\^y (75) 

Then a flux larger than 1 particle in T years can be achieved if 

(^^y"^U^y\AT>10^ (76) 

\lTeVJ \lkpc) ^ ' 

where A is the detector area in m^ . A detector of effective area lO^m^ (compare with 400m^ for 
H.E.S.S.) will therefore receive more than one particle/year if m w O.lTeV and D w Ikpc and t^l. 
Obviously, both e and the fraction of particles effectively detected can be much less than unity, 
especially if the particles are weakly interacting; the quick estimate of this Appendix is meant only 
to give an idea of the orders of magnitudes. 
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Figure 7: Proper motion geometry. The source moves with velocity Vs from S' to S. When the source is at 
S' it emits a particle with velocity v; when is in S it emits light or a ultra-relativistic particle with velocity 
c. The distance OC is the impact radius ^pm. The time it takes for the particle to travel from S' to the 
observer must equal the times it takes for the source to move from S' to S and for the relativistic signal to 
propagate from S to the observer. 



Appendix C. Source's proper motion 



If the source has a non-negligible proper motion with respect to us then the non-relativistic 
particles that we observe today were emitted from a position that depends on their velocity. The 
particles will have therefore a different arrival direction with respect to the relativistic counterpart. 
This will spread the particles on the plane that contains the source velocity vector and the observer, 
inducing a "proper motion alignment". 

Let us now denote with 6c (see Fig. [7]) the angle between the source velocity v^ and the direction 
of arrival of a relativistic signal (photon or particle) and with 9 the corresponding angle for a 
non-relativistic particle and let us assume a uniform motion. All velocities are in the observer rest 
frame. In order for the two signals to arrive at approximately the same time we require 



1 



cot( 



1 



COt^r 



V smi 



(77) 



for non-relativistic Vg = |v<;|. This condition can be solved analytically to obtain 9{0c,Vs,v). It 
is more useful however to derive an approximate solution. A non-relativistic particle will arrive 
therefore with a proper motion deviation angle 



A 



PM 



with respect to the relativistic signal. For small deviations Apjv/ we obtain 

{l — v)vsSm9c VsSm9c 



^PM 



w, cos 9r 



V ~ Vs cos 9c 



(78) 



(79) 



where the last step applies for 1; <C 1. Note that v is always larger that Vs cos 9c due to the 
composition of velocities, so ApM > 0. Finally, if the source is moving slowly with respect to the 
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particle velocity (a reasonable assumption since the velocity dispersion in the Milky Way at the 
Sun's location is around 50-100 km/sec (see e.g. |13|) we obtain 

APM « !^^ . ^ (80) 

V Ts CV 

in terms of the impact parameter ^pm and the distance to the source rg or in terms of the proper 
motion v = cVsSiii(9c) /rg (the velocity of light appears explicitly now because v is the angular 
velocity). Unless the source proper motion is negligible, the particles will arrive to the observer 
aligned on a direction which is intermediate between the lensing direction (source-lens) and the 
source proper motion vector. Assuming a proper motion of /xarcsec/year we have 

ApM«(9deg),(^)(^)- (81) 

i.e. a typical deviation of ten degrees from the direction induced by the non-relativistic lens. In 
general the resulting detection trail will be a combination of the lensing and proper motion effect, it 
will point in a random direction and for large v it will no longer be an exact straight line due to the 
different behavior with respect to u. In principle however the parameters /i, r^ of a given galactic 
source can be estimated through direct astrometric observation and utilized to subtract the proper 
motion from the detection map. 

Let us finally notice that the particle velocity v and arrival angle in the earth rest frame are 
different from the corresponding quantities in the source rest frame. Suppose the source has a 
peculiar velocity vector Vg with respect to us. A particle is emitted with velocity u in the frame 
in which the source is at rest, and received with velocity v in our rest frame. All velocities are in 
units of c. The relativistic addition theorem for velocities says that 

v = 7-^ + -fi- (82) 

1 — im 7s (1 — vu) 



where 7^ = (1 — v1) ^1"^ and 



Putting 



we find 



U|| = Vs(vs-u)/y2 (83) 

uj_ = u - U|| (84) 



u = u{cos6l',sin6l'} (85) 

V = u{cos6', sin6'} (86) 



. ^ usa\Q' ,„„. 

sin 6* = ^ (87) 



h^(vs +ucos0'y + u^ sin^ 9')] 
7s (1 -I- uvg cos 9') 

The equation for 9' can also be written as 

9 usin9' 

2 "fs{vs+ucos9') + [j^{vs + ucos9')^ + u'^sin'^9')]^/^ 

which for u — >■ 1 becomes the standard aberration equation 

« 9'/l-z;,^^/^ 



(and at the same time u — v). 



tan-.tan-^^) (90) 



